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The data on pp elastic scattering at 1.8 and 1.96 TeV are analysed in terms of real 
and imaginary amplitudes, in a treatment of very high accuracy, covering the whole 
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slopes. A method is introduced for determination of the total cross section and the 
other forward scattering parameters and to check compatibility of E-710, CDF and 
the recent DO data. Slopes Br and Bj of the real and imaginary amplitudes, treated 



'-£3 | as independent quantities, influence the amplitudes in the whole t-range and are 

O ' important for the determination of the total cross section. New necessary generalized 

D. 

_C ■ expression is derived for the Coulomb phase. Comparison is made of the values of 
a and Bj obtained in forward and full-t approaches, with a new determination of 
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qq \ these quantities. The amplitudes are fully constructed, and a prediction is made of 
00- 

"^J- ■ a marked dip in da/dt in the \t\ range 3-5 GeV due to the universal contribution 



of the process of three gluon exchange. 



INTRODUCTION 



The precise knowledge of total cross section and scattering amplitudes in pp and pp 
elastic scattering at high energies is essential for understanding the QCD interactions and 
hadronic structure, and also for the parametrization and extrapolation of the total cross 
section that may pass through the LHC measurements and go up to the study of ultra-high 
energy phenomena in cosmic rays rays [lj. However, at high energies, due to the smaliness 
of the ratio p of the real to the imaginary parts of amplitudes at t = 0, together with the 
absence of data for small \t\, turn the extrapolations toward to \t\ — > very delicate. It is of 
fundamental importance to characterize well the \t\ — > limits of the scattering amplitudes 
that are used to determine forward slopes and total cross section. 

It is universally understood that the real and imaginary amplitudes in pp and pp elastic 
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scattering amplitudes reflect the non-perturbative QCD dynamics, determined by overall 
features of the proton and antiproton structures. Regge-like behavior characterizes the s 
and t dependences at large s and small \t\. As a consequence, there appears a dip or an 
inflection point in the differential cross section da/dt, due mainly to a zero in the imaginary 
part, and the detailed shape around this region is influenced by the magnitude, sign and 
form of the real part. An analysis of the interplay of real and imaginary amplitudes is 
necessary to explain this behavior [2J. 

In this work, we present a high precision description of the data on pp elastic scattering 
at 1.8 and 1.96 TeV, consistently extended from the forward to the backward regions. This 
is particularly opportune in view of the publication of the new measurements at 1.96 TeV by 
the DO Fermilab experiment sj]. We introduce analytical forms for the amplitudes that are 
not pure exponentials, and identify their zeros, signs, ranges of dominance and the interplays 
that fix the details of the observed inflection, dips and bumps. The analytical form of our 
amplitudes are based on previous studies made in the framework of the Stochastic Vacuum 
Model k|. 

The present work applies previous studies, enriched with the controls offered by dispersion 
relations for slopes [5] and with the incorporation of new data. Here, we couple the previous 
1.8 TeV data with the recent measurements at 1.96 TeV and verify that they are compatible. 

It is stressed here that in the analysis of the data, we use information from forward 
dispersion relations to control parameters of the full description and particulary emphasize 
the importance of the difference of the slope paramters, Br and Bi, of the real and imaginary 
parts. It is commonly assumed in the analysis of data that these slopes are the same, but 
this is wrong theoretically, and this fact is often overlooked due to the smallness of the p 
parameter. However, when we want to describe consistently the scattering amplitude for 
the full \t\ range, this difference is crucial. Such a description which covers the long t region 
constrains the quantities of the forward range. This is particularly true and important for 
the real amplitude that is small with respect to the imaginary part in the forward direction, 
but not at large — t — . Actually, the differential cross section at high \t\ is dominated by 
the real part. 

On the other hand, in the very large \t\ domain, the perturbative QCD effects become 
dominant. We investigate the connection of the measured points at 1.8-1.96 TeV with data 
of the perturbative tail of large \t\ measured at 27 GeV [6]. It is known that this tail is 
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energy independent and formed by a real amplitude due to three-gluon exchange |7|, with 
opposite signs for pp and pp (positive in pp and negative in pp) scattering. We show that 
inclusion of such a universal perturbative amplitude leads to a striking prediction for the 
behavior of the cross section. In pp scattering, when added to the non-perturbative positive 
real part, the perturbative term creates a third zero, located in the region about 3-4 GeV 2 . 
As the imaginary part is less important in this domain, a marked dip is formed with this 
cancellation. 

Our treatment of the whole data with overall high precision leads to definite prediction 
of the forward scattering parameters in the context of one analytical form. We thus have 
a determination of the total cross section a and of the quantities p, Bj and Br (see Eq. 
( IT5l) ) . This treatment is compared to the analysis of a subset of forward points in terms 
of purely exponential amplitudes, in which we obtain different values for a and Bj (see Eq. 
(|T6|) ). We show that the two determinations are statistically equivalent, and call attention 
that the discrepancy of the two approaches would be reduced if more data at low — t — were 
available. Thus the determination of o remains model dependent. Obviously, we avocate 
that the determination based on the full data is more interesting, since it incorporates 
properties of curvatures and zeros of the amplitudes. 

We organize the present work as follows. In Sec. 2, we present scattering amplitudes that 
describe the whole |£| range , and the necessary quantities are defined, with a discussion of the 
role of the universal perturbative amplitude for large \t\. In Sec. 3, the precise analysis of the 
data at 1.8-1.96 TeV to determine all quantities is presented. In Sec. 4, we summarize our 
results and discuss the prediction introduced in this work. In the Appendix, we derive new 
expression for the Coulomb phase adapted to incorporate the difference in slope parameters 
of real and imaginary amplitudes. 

II. GENERAL DESCRIPTION OF FORWARD AND FULL t SCATTERING 

In the treatment of elastic pp and pp scattering in the forward direction, with amplitudes 
approximated by pure exponential forms, the differential cross section is written 
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where t = —\t\ and we have assumed different vaiues for the slopes Bj and Br of the 
imaginary and real amplitudes. In the following discussion, we use the unit system where a 
is in mb(milibarns) and energy in GeV, so that (hey = 0.389 mbGeV 2 . 
The Coulomb amplitude F c (s,t) enters for pp/pp with the form 

F c (s, t)e- •<*.*) = (-/+) •(••*) F p 2 roton (t) , (2) 

where a is the fine-structure constant, $(s,£) is the Coulomb phase and the proton form 
factor is written 

F piotoa (t) = [0.71/(0.71 + |t|)] 2 . (3) 

In usual treatments of the data, no distinction is made between Br and Bj slopes, in 
contradiction with expectations from dispersion relation |5J, and a necessary generalization 
of the commonly used assumption is required. However, when Br ^ Bj, treatment of the 
Coulomb interference requires a more general expression for the Coulomb phase, which is 
derived in the appendix. 

In elastic pp and pp scattering at all energies above y/s = 19 GeV, the real and imaginary 
amplitudes have zeros located in ranges \t\ ~ (0.1 — 0.3) GeV 2 and \t\ = (0.5 — 1.5) GeV 2 
respectively, and the use of exponential forms beyond a limited forward range leads to 
inaccurate determination of the characteristic forward scattering parameters a, p , Bj and 
Br. To obtain precise description of the elastic da/ dt data for all \t\ , we introduce amplitudes 
with forms I2f| 



T R (s,t) = a R (s)exp{-l3 R {s)\t\) + X R (s)^ R (j R (s),t) + ^F c (t) cos(a$) 



and 



Tj(s,t) = aj(s)exp(-/?j(s)|*|) + Aj(s)ttj(7/(s), t) + V^F c (t) sin(a$) 
with the shape functions 



#jc(7*(s),*) = 2e 
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(4) 
(5) 

(6) 



a/1 + a |£| a/4 + a \t\ 

where K = R for the real amplitude and K = I for the imaginary amplitude. We here 
have quantities «/, (3i, 77, \j, a R , (3 R , j R , X R whose roles are explained below. 7^ is 
dimensionless, while ax, Jk and (3k are in GeV -2 . These forms have been developed in 
the application of the Stochastic Vacuum Model to pp and pp elastic scattering [2| , and the 
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fixed quantity a® = 1.39 GeV -2 is related to the square of the correlation length of the gluon 
vacuum expectation value ( a = (0.2 — 0.3) fm) |4|. 

From the above expression, we can express the total cross section a (s) ,the ratio p of the 
real to imaginary amplitudes, the slopes B R j of the amplitudes at t = 0, and the differential 
cross section da/dt as 

a(s) = A^(hcf (aj(s) + Aj(s)) , (7) 

= T R (s,t = 0) = a R (s) + X R (s) 
P[S) ~ Tj(s,t = 0) ~ ai (s) + \i{s) ' 1 ' 



B K (s) 



7V(s, £) dt 
1 



1 dTx(s,t) 

a K {s)(3 K {s) + -X K (s)a (6j K (s) + 7)] , (9) 



t=o 



oc K (s) + A A '0 
da 



^ = (M \T R (s,t)+iTj(s,t)\ 2 . (10) 

We have thus defined the form of the amplitudes for all t at each energy. The parameters 
must be obtained from phenomenological analysis of the data, with control from dispersion 
relations for amplitudes and for slopes. The forms of Eqs. (jlj), @ ar e able to describe the 
imaginary and real amplitudes at all energies, with consistency in their features (magnitudes, 
signs, locations of zeros), and with smoothness in the energy dependence of the parameters. 
Values of p and B R must be related with a and Bj respecting dispersion relations. The 
forms have been used successfully to represent the data at all energies 2]. 

This description represents the nonperturbative QCD dynamics that is responsible for 
soft elastic hadronic scattering. They effectively account for the terms of Regge and eikonal 
phenomenology that determine the process up to \t\ ranges up to about \t\ ~ 2.0 GeV 2 . 
These forms disentangle the imaginary and real amplitudes and are able to reproduce with 
good accuracy all da/dt behavior. 

In the applications of these forms to phenomenology it was found that the imaginary 
amplitude presents one zero located in the range (0.5 - 1.5 GeV 2 ). The real amplitude 
presents one zero at low \t\ (\t\ < 0.3 GeV 2 ) , according to a theorem by Martin B|, and 
a second zero whose location determines the shape of da/dt around the dip (or inflection 
point). As a general behavior, we have that the imaginary part Eq. (jHJ) is negative and the 
real part Eq. (j3J) is positive for \t\ larger than 1.5 GeV 2 . These features are simple and all 
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data are described accurately. The amplitudes are draw in the LHS of Fig [T] , with marks 
at their zeros. 

It is observed that after the dip (or inflection point) the behavior of the differential cross 
sections becomes increasingly energy independent. The elastic pp experiment at \/s = 27 
GeV [6j has measured the range from 5.5 to 14.2 GeV 2 and is the only experiment covering 
such high values of — t — . This distribution at high \t\ shows remarkable universality: 
at all energies y/s = 23.5, 30.7, 44.7, 52.8 and 62.5 GeV , namely at all energies where 
measurements have reached the intermediate — t — region, da/dt approaches the same set 
of points of the 27.4 GeV experiment. 

The observed da/dt at the tail has a dependence of form l/|t| 8 , and has been explained 
by Donnachie and Landshoff 7j] as being of perturbative origin, due to a contribution of 
three gluon exchange. This term is real and has an amplitude of the form 
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m 2 (\t\) 



where 

«.(!<!) = — t / ( 12 ) 



1 / m 2 (\t\)+t\ 

log ( J 



is the strong coupling constant and ?7t(|t|) is the gluon effective mass {{J. The factor 3 in 
t ~ (y^/3) 2 comes from the assumption that each gluon carries one third of the momentum. 
The normalization factor N is negative and determined by the nucleon structure. To extend 
our description to include the very high \t\ range of this form, we introduce a term R ggg (t) 
, writing 

TVaii)(M) = a R (s)ex P (-^ R (s)\t\) + X R (s)^ R ( lR (s),t) + V^F c \t)cos(a^)+R ggg (t) , (13) 

where R ggg (t) is chosen so that the differential cross section to be dominated by a term of the 
form |i|~ 8 for large \t\ values (say above 2.5 GeV 2 ), while for small \t\ the amplitude stays 
unaltered as the original nonpertubative expression. The perturbative three gluon exchange 
has opposite signs for pp and pp scattering, being positive for pp and negative for pp. We 
then take the following expression, 

R ggg {t) = ± 0.45r 4 (l-e- a005 l'l 4 )(l-e- ai l*l 2 ) , (14) 

where the signs ± apply to the pp and pp amplitudes respectively. The factor 0.45 is chosen 
to reproduce the Faissler measurements and the last two factors are written to suppress 
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FIG. 1: Amplitudes in pp and pp elastic scattering shown in different ranges and scales, described 
by Eqs. ([!]), ©, (J6]) and determined by phenomenology. In the |t| range up to about 2 GeV 2 the 
amplitudes are governed by nonperturbative dynamics and are qualitativelly similar for the pp and 
pp, with one zero for Tj and two zeros for Tr. Tj remains negative and goes fast to zero, while 
at \t\ ~ 3 GeV 2 the nonperturbative Tr is positive and dominates. In pp scattering the negative 
contribution of the 3 gluon exchange term cancels Tr, forming a third zero and a marked dip in 
da/dt. 

smoothly the perturbative contribution for small \t\. The transition region from 2 to 5 
GeV 2 contains information on the nature and interplay of nonperturbative and perturbative 
contributions, which must be investigated. 

In the analysis of the scattering amplitudes in the nonperturbative domain, we find that 
the real amplitude is positive in the transition region, for both pp and pp systems. In pp 
scattering, the inclusion of the negative tail amplitude, cancels the positive nonperturbative 
part, creating a third zero. Then, as the imaginary part is not dominant, a marked dip may 
be observed in da/dt in the transition region. This is shown in Fig. [31 
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III. ANALYSIS OF ELASTIC pp DATA AT 1800 GEV 



The available data at 1800 GeV consist of 51 points from the E-710 experiment 10| for 



0.0339 < \t\ < 0.627 GeV 2 and 26 points from the CDF experiment for 0.26 < \t\ < 1.26 
GeV 2 . Recent data at 1.96 TeV [3] from the DO experiment with 17 points with 0.26 < 
\t\ < 1.26 are considered to supersede part of the E-710 data , forming a set of 35+17=52 
points, that we call Set I. The new data from the DO experiment at 1.96 TeV have here 
been reduced by a factor (1.8/1.96) 3232 = 0.973 based on Regge phenomenology in order to 
adjust all data to the same 1.8 TeV energy ( as the t range involved is small we neglect the 
\t\ dependence of this factor). We merge with these the 26 points from the CDF experiment, 
that covers about the same range of the 17 DO points at 1.96 TeV, mounting Set II with 78 
points. 

The existing data stops at about \t\ = 1.3 GeV 2 , leaving the higher \t\ region without 
information and without constraints to fix the connection with the range of the perturbative 
tail. Thus the parameter (3r cannot be fixed uniquely. As its value is crucial for the 
prediction of the position and depth of a dip in the transition region for 1.8 TeV, we present 
two alternative choices, with /3r = 1.15 and 1.40. The two solutions have the same values for 
aj, (3j, Xj, (and thus a), p, £>/, B R . Noting that 77-, 7^, a R are fixed from other quantities, 
the fitting of each solution is needed only to provide A^. The results are given in Table [J 
where we observe remarkable quality in the description of the data, with very low < x 2 > 
values, smaller than 1. 

In a parametrization covering all \t\ values using Eqs. (HI), (jSJ), (EJ), we find that sets I 
and II are compatible with each other, leading to approximately the same numerical values 
of the parameters. The < x 2 > values are different in the two cases, because of redundancy 
of equivalent points (the 26 CDF points and the 17 DO points are show to be statistically 
equivalent in their central values) and because of the larger error bars of the CDF data. Thus 
the more complete Set II has an improved average < x 2 > = 0.67 compared to < x 2 > = 
0.80 of Set I. 

In what follows we adopt the parameter values obtained with the more complete Set II, 
with N = 78 points. The values of the physical quantities are 

a = 72.75 ± 0.18 , p = 0.138 ± 0.004 , B 1 = 16.80 ± 0.04 , B R = 26.23 ± 0.37 . (15) 

We recall that the above analysis is based on our analytical expressions applied to all 
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TABLE I: Quantities of the all-t representation. Common values : a = 72.75 ±0.18 , p = 
0.138±0.0004 , Bi = 16.80±0.04, B R = 26.23±0.37 , m = 12.9632±0.0248, A/ = 13.4156±0.0399, 
Pi = 4.2060 ± 0.0082. Set I built with E-710 and DO data. Set II has additionally the CDF data. 
The real amplitude parameters /3r and Xr characterize different solutions. < x 2 > is the average 
of squared relative theoretical/experimental deviations. 



Set 




Pr = 1-15 




Pr = 1-40 


I 


Xr 


= 3.5714 ± 0.0094 


Xr 


= 3.5628 ±0.0105 


N=52 




X 2 = 0.80 




X 2 = 0.80 


II 


Xr 


= 3.5628 ± 0.0088 


Xr 


= 3.5530 ± 0.0099 


N=78 




X 2 = 0.67 




X 2 = 0.67 



\t\ . We believe that this analysis leads to a reliable determination of the so called forward 
scattering parameters a, p , Bj , Br that are characteristic of forward scattering. Our 
analytical forms intrinsically have curvatures in the amplitudes that allow them to cross 
zero so that they are able to describe the characteristic structure of da/dt, with dips, bumps 
and inflection points. 

Now, in order to check the sensitivity of our fits of above parameters to the global form 
of our amplitude, we analyse separately the forward scattering region that usually serves 
to the determination of a, p, Bj. In such analyses, traditionally the amplitudes at low \t\ 
are described with pure exponential forms as in Eq.([T]), and must be restricted to a low \t\ 
range. 

The 1.8/1.96 TeV data are poor in the forward direction, and do not include the interval 
from 10" 3 to lCr 2 GeV 2 , with the E-710 data starting at \t\ = 0.0339 GeV 2 . The CDF data 
are restricted to the intermediate |£| range, \t\ ^ 0.26 GeV 2 , and do not give good basis for 
extrapolation to find the optical point and determine the total cross section. 

We select the 24 first points 0.0339 < \t\ < 0.0827 of the E-710 experiment and fit with 
4 parameters using Eq. ([T]). Using Eq. Qj, our results, with < x 2 >— 0.62, are 

a = 71.92 ± 0.23 , p = 0.138 ± 0.031 , Bi = 16.20 ± 0.12 , B R = 26.23 ± 9.27 . (16) 

Thus the parameters of the imaginary amplitude (a and B{) are not the same as in the 
long \t\ analysis in Eq. (ll5D and the larger error bars in p and Br must be noted. The lines 
representing the two kinds of description of the data are shown in the RHS of Fig. [2j 
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FIG. 2: LHS: Data in pp scattering at 1.8 TeV : 35 points from A.Amos et al 
points from V.M. Abazov et a 
26 points from F.Abe et al 



with 17 



(slightly converted from 1.96 to 1.8 TeV) in open circles, and 
111 ] in black circles. In the LHS solid (and another superposed line) 
show analytical forms that describe the data (see Table H]). RHS: the first 24 forward points, with 
\t\ < 0.1 GeV 2 treated with forward scattering form of Eq. ([1]) in dashed line; solid line represents 
the same full-t line of the LHS. There is no visible distinction in the range of the experimental 
points, but the limits \t\ —> lead to different values of a and Bj, given in Eqs. (I15p and (|16p . 



We can compare numerically the quality of the two descriptions of the 24 data points. 
The average squares of relative deviations of theoretical values and data are < %2 >= 0.606 
and < x2 >= 0.613 for the forward and the long-i forms, respectively. We thus see that the 
two representations (forward and long-t forms) of these data set are practically equivalent 
in this domain. Different values of a (72.75 mb and 71.92 mb) and of Bj (16.80 GeV~ 2 
and 16.20 GeV -2 ) indicate the effect of the use of two different approaches to the idealized 
\t\ = limit. 

We must notice that the real amplitude parameters p and Br are the same (values 
suggested by dispersion relations) in Eqs. f lT5|) and ffTB"]) . because they do not receive strong 
influence from the points of the set of 24 forward points. Actually these parameters become 
important in the determination of the properties of the real amplitude only for larger \t\. 

The ambiguity in the determination of a and Bj (compare Eqs. ( llSp and ( fTBl ) occurs 
because the data are poor, inexistent below 0.0339 GeV 2 . It is clear for us that the value 
of a reported as output (around 80 mb) from the analysis of the CDF experiment based 



11 



TABLE II: The inflection point in da/dt, positions of zeros of the amplitudes, locations of the 
predicted dip and bump, and ratio characterizing the shape of this dip. All quantities in GeV 2 , 
except the dimensionless ratio. Integrated elastic cross section a e \ = 16.65 mb 



GeV" 2 


|* inn 

GeV 2 


{da/dt)irft 
mb/GeV 2 


Zi Z R (1) 
GeV 2 GeV 2 


Zr(2) Zr(3) \t dip t bump ratio 
GeV 2 GeV 2 GeV 2 GeV 2 


1.15 


0.7370 


0.010 


0.6240 0.1774 


1.4447 3.6955 3.7376 4.6869 9.0413 


1.40 


0.7356 


0.010 


0.6240 0.1779 


1.5744 3.0633 3.3198 4.1457 2.0319 



on points \t\ > 0.26 GeV 2 must be revised. The energy 1.8/1.96 TeV is important for the 
parametrization of a(s) and its extrapolation for ultra-high energies. We note that the 
value derived from the CDF data without desentangling the real and imaginary amplitudes 
gives an increase in the a(s) curve, whereas our value enters more smoothly in the global 
behavior of a(s). Obviously we recommend the use of a long \t\ description, with proper 
analytical expression, that can be different from ours, but must have the similar features for 
the amplitudes such as proper magnitudes, curvatures, zeros and signs. 

An additional remark about this question of analytical shapes and roles of limits comes 
from the observation that if we put the values for a = 72.75 mb, Bj = 16.80 GeV -2 (together 
with p and Br values) in the usual forward amplitudes forms, Eq.([T]) , we obtain x 2 — 0.744 
for the 24 points, that is formally worse than 0.606 of the direct fitting. We thus read the 
concept: in a case like this, of poor data at low \t\, the parameter values represent \t\ — > 
limits of specific analytical forms, and are to be used associated to the expressions for which 
they have been found. Thus the values of a and Bj may be considered as model-dependent. 
Eq. (Q]) is not a better model. Actually, with Br = Bj as usually made, it is wrong. 

We recall values of the scattering parameters that are found in papers by experimentalists 



12|: a = 71.71 ± 2.02 mb , a el = 15.79 ± 0.87 mb, B = 16.3 ± 0.3 GeV" 2 (and also 
B = 16.98 ± 0.22 GeV" 2 ) , p = 0.135 ± 0.044). We must remark that B (slope of da/dt) is 
not the same as Bj , and that Br ^ Bj gives a different meaning to the parameter p. 

In Table [Til we give the values of \t\ at the zeros of the amplitudes, and the locations 
of the dip and bump in da/dt due to the contribution of the 3-gluon exchange term. The 
quantity ratio = ( da/dt) bump /(d<r/dt)di P that informs about the shape of the dip depends 
strongly on the values of the parameter (3 R , that must be determined by experiment. 
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IV. SUMMARY AND DISCUSSION 

In the present work, we have presented a very precise description of the elastic scattering 
amplitude for the pp collisions, merging the recent 1.96 TeV and the former 1.8 TeV data. 
We introduce analytical forms for the amplitudes in the full range of \t\, and identify their 
zeros, signs, ranges of dominance and the interplays that fix the observed details. We believe 
that with our method a more correct determination of the forward scattering parameters at 
1.8 TeV can be obtained, together with the global form of the amplitude shown in Fig. [TJ 

We stress that to determine the total cross section and related observables, we must 
study the properties of the amplitudes, and not a simple phenomenological presentation of 
the squared observables as da/dt which may be misleading : in our practice we found cases 
of very low \ 2 values for the description of da/dt, obtained from unacceptable amplitudes. 

Figj2] clearly shows how delicate is the extrapolation of experimental data of da/dt to- 
wards \t\ = 0. Our work revises the values of total cross section and slope parameters that 
are reported in the literature, suggesting new values. 

Our analysis indicates that at 1.8 TeV the imaginary zero is in a region where the real 




FIG. 3: Predictions for the contributions of real da R /dt and imaginary da 1 /dt parts of da/dt in 
the presence of the real perturbative tail due to 3 gluon exchange. In pp scattering the negative 
sign of the tail causes a zero in da R /dt and a dip in da/dt located in the range 3-5 GeV 2 . The 
RHS figure shows two examples of the dip structure, formed with /3r = 1.15 (solid) and /3r = 1.40 
(dashed) as given in Tabled 
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part is not weak. Thus there is an inflection point, but not a marked dip as occurs in pp at 
lower energies. With the increased energy, the imaginary zero has been displaced to lower 
— t — values more rapidly, and is now distant from the second real zero. The real part is 
positive and much smaller at t = 0, and passes through zero at \t\ ~ 0.2 GeV 2 , becoming 
negative until it crosses zero again, after the inflection point that occurs near 0.74 GeV 2 . 
The second real zero occurs at \t\ wl.5 GeV 2 , while the imaginary part remains always 
negative after its zero near 0.6 GeV 2 . 

For large — t — this real part becomes increasingly dominant over the (negative) imag- 
inary part, that has a rapidly decreasing magnitude. This is the range where the real 
nonperturbative amplitude dominates. This behavior is observed in Fig. (TJ with the line Tr 
with no tail. We are here in the region of transition between non-perturbative and pertur- 
bative dominances. As the nonperturbative real part vanishes, the tail remains, giving to 
the differential cross section the characteristic shape l/|t| 8 - This behavior of dip formation 
is predicted equally for the pp cross sections at 541 and 1800 GeV. 

The situation is different for pp at 52.6 GeV, where the second zero would occur at a 
higher \t\ (about 2 GeV 2 ) and is not formed due to the action of the negative perturbative 
tail. At 541 and 1800 GeV the real amplitude become positive before the negative tail enters. 

To conclude, we propose an experimental investigation of da/dt in elastic pp scattering 
in the collected data of the DO collaboration at values of \t\ beyond those already published, 
where a marked dip may appear. Typical expected results are show in Fig. [3] An observed 
structure will identify the sign and magnitude of the real amplitude in this region. It becomes 
clear that the investigation of the transition range from nonperturbative dynamics to the 
tail region dominated by perturbative contributions is very important. 



V. APPENDIX: THE COULOMB PHASE 



Here we derive an expression for the phase appropriate for cases with Br ^ Bj. 
The starting point is the expression for the phase obtained by West and Yennie 

r0 dt' r. F N ( s ,try 



= (-/+) 



ln( -~] 



ip 2 



w -t 



1 



F N (s,t) 



(17) 



where the signs (— /+) are applied to the choices pp/pp respectively. The quantity p is the 



proton momentum in center of mass system, and at high energies 4p 2 



s . 
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For small \t\, assuming that F N (s,t') keeps the same form for large \t'\ (this approximation 
should not have practical importance for the results), we have 

F N (s, t') Fg(s, 0)e B «*'/ 2 + % Fj N (s, 0)e B ' t '/ 2 



F N (s,t) ' Fg(s,0)e B ^ 2 + t Fj N (s,0)e B ^ 2 



c + i 



,B R (f-t)/2 



+ 



3 Bl(t'-i)/2 



c + i 



where 



pe 



{B R -Bi)t/2 



(18) 



(19) 



The calculation is explained in detail in the appendix. The integrals that appear in the 



evaluation of Eq. ffTTl) are reduced to the form 

"° _dt' 

-Ap- 



1(B) 



1 _ e B{t'-t)/2 



(20) 



that is solved in terms of exponential integrals 15] as 



— + In — 

2 J L 2 



(4p 2 + t)]+ln[-^]+2 7 . (21) 



The real and imaginary parts of the phase are then written 



= ("/+) 



In 



.s/ + c 2 + 1 



c 2 I(B R ) + I(Bj) 



and 



= (-/+) 



c 2 + l 



I(Bj) - I(B R ) 



(22) 



(23) 



With a in mb and t in GeV 2 , the practical expression for da/dt in terms of the parameters 
a, p , Bi and B R is 

,B R t/2 n 2 



da 

— = it (nc 
dt v 



p a e 



+ F G e Q * / cos(a$ i? ) 



4vr (he) 2 

At high energies and small \t\ we simplify 



+ 



a e 



47T (ftc 



Bjt/2 -, 2' 

+ i^e"* 1 sin(a$ fl ) 



• (24) 



4p 2 + t -)> s 

and then the functional form of /(-B) is written 



/(BH£l( £f)_ £ , ( _ f)+lll(f)+ln( 



Bt\ 

-2) +2 ~« 



(25) 
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FIG. 4: Comparison of the Coulomb phase a&R from our calculation of Eq. (|22|) with the West- 
Yennie WY expression in Eq. (I17p . The unrealistic large values of parameters p = 0.3 and 
Br/Bj = 5 are chosen in order to enhance differences. The calculation is made at y / i=50 GeV, 
but the values of a&n(s,t) do not show explicit dependence on s up to the LHC energy. The plot 
of a$/ in the RHS shows extremely small values (notice the scale), about 100 times smaller than 
a&R , so that we can safely put a<3?/ equal to zero in Eq. ()24jl . 
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The expression for the phase can be compared with the usual expression from West 
and Yennie 



WY 



("/ + )« 



7 + In ( 



Bt\ 



(26) 



A numerical comparison is presented in Fig. HI where we choose strong unrealistic parameter 
values p = 0.3 , B R /Bj = 5 to enhance possible differences, and the lines are drawn for 
Bj = 15 GeV -2 , with the da/dt slope B of the WY formula taken as given by 

P 2 B R + Bj 



B 



(27) 



1 + p 2 ' 

At the high energies of our interest, the phase is insensitive to s. 

The plot of a$/(s,t) in the RHS of Fig. H] shows that its values are extremely small so 
the the imaginary part of the phase can be safely put equal to zero. 
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